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ABSTRACT 

We extend investigations of mass entrainment by jets, which previously have 
focused on cylindrical supermagnetosonic jets and expanding trans-Alfvenic jets, to a 
set of expanding supermagnetosonic jets. We precess these jets at the origin to excite 
the helical mode of the Kelvin-Helmholtz (or KH) instability, in order to compare 
the results with predictions from linear stability analysis. We analyze this simulation 
set for the spatial development of magnetized mass, which we interpret as jet plus 
entrained, initially unmagnetized external mass. As with the previous simulation sets, 
we find that the growth of magnetized mass is associated with the growth of the KH 
instability through linear, nonlinear, and saturated stages and with the expansion of 
magnetized material in simulated observations of the jet. From comparison of measured 
wavelengths and wave speeds with the predictions from linear stability analysis, we 
see evidence that the KH instability is the primary cause for mass entrainment in 
these simulations, and that the expansion reduces the rate of mass entrainment. 
This reduced rate can be observed as a somewhat greater distance between the two 
transition points separating the three stages of expansion. 



Subject headings: galaxies:jets — instabilities — MHD 



1. INTRODUCTION 

Mass entrainment of an external medium into a jet flow affects the spatial development of 
the jet and may be a crucial factor in the division between FR I and FR II (Fanaroff & Riley 
1974) radio sources. For example, deceleration, presumably via mass entrainment, of a turbulent, 
supersonic jet is one of the basic assumptions in the model of Bicknell (1994, 1996) that reproduces 
the break between FR I and FR II sources in the radio-optical plane. Additionally, there is 
evidence for jet slowing in the low-power radio source B2 1144+35, which is observed to decelerate 
from 0.95c to 0.02c between 20 pc and 24 kpc from the core (Giovannini et al. 1999). Mass 
entrainment may be important even in slower jet flows, such as the much slower (0.1c) parsec-scale 
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jets in Seyfert galaxies (Ulvestad et al. 1999). Since our simulations do not include radiative 
cooling, which is a crucial element in the energetics of protostellar jets, the simulations we present 
here are not strictly applicable to those systems, although perhaps they are appropriate for the 
jets in Galactic super luminals. 

Rosen et al. (1999, hereafter RHCJ) confirmed that a reasonably strong magnetic field can 
slow mass entrainment by a three-dimensional, nonrelativistic, cylindrical supermagnetosonic 
jet. RHCJ showed that the spatial growth of mass entrainment is a three-stage process that 
corresponds to the growth of the Kelvin-Helmholtz (KH) instability, from a linear stage through a 
nonlinear stage to a stage where the growth has become saturated. RHCJ also found three stages 
of growth in the width of simulated radio intensity maps of the jets, and that these three stages in 
the width are associated with those of the growth of the jet plus entrained mass. However, some 
detailed differences did exist in the position of the transition between stages as a function of jet 
density. These differences usually occurred for jets lighter than the external medium, and took the 
form of continued mass entrainment even though the width in a simulated intensity image of the 
jet remained constant. A jet that entrains mass could have an observable spine-sheath morphology, 
in which the spine is more disrupted in a jet that is lighter than the external medium. In fact, 
recent application of the "spectral tomography" process (Katz-Stone et al. 1999) has revealed the 
presence of a low intensity, steep spectrum sheath surrounding a high intensity, flat spectrum core 
or jet spine in a few wide-angle tail radio sources. In addition, RHCJ found that the linear growth 
of jet plus entrained mass was related to a linear progression as higher order, small amplitude KH 
modes with short growth lengths are overwhelmed by lower order, larger amplitude KH modes 
with larger growth lengths. 

In this paper, we extend our mass entrainment analysis to conical supermagnetosonic jets 
that have inlet conditions similar to some of the jets in RHCJ. Based on the KH linear stability 
analysis for supermagnetosonic jets in Hardee, Clarke, & Rosen (1997, hereafter HCR), the growth 
length at maximum growth (i.e., the inverse of the spatial growth rate at maximum growth) is 
proportional to the radius of the jet. Additionally, the maximum growth rates are reduced along 
an adiabatically cooling jet, in which the Mach number increases with the jet expansion (Hardee 
1986). Thus, in general on supermagnetosonic jets the instability should grow more slowly as the 
radius expands than on a cylindrical jet. We investigate whether expanding supermagnetosonic 
jets are more stable than similar cylindrical jets and also whether expanding jets show the similar 
association between the spatial development of mass entrainment, the spatial development of 
the KH instability and the spatial development of the apparent width of magnetized material 
in simulated radio intensity images that has been seen previously in cylindrical jet simulations 
(RHCJ). 

Examination of the potentially smaller effects of the KH instability on conical jets is an 
important extension of previous work, as there are observational examples of conical jets, some 
of which expand by orders of magnitude and remain undisrupted for more than 100 kpc. For 
example, the jet in M87 expands at a nearly constant rate in a recent high resolution (O'/l « 8.2 
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pc) radio image of the source out to knot A (Biretta, Zhou, & Owen 1995). Another example of 
jet expansion is provided by NGC 6251, which has a jet that passes through a series of expansions, 
each of which has an average half-opening angle of ~4° for at least 10 kiloparsecs in length (Perley 
et al. 1984). The jets in both of these radio sources appear more laminar than turbulent in the 
regions where the jets are conical. 

Since we aim to investigate the spatial growth of the KH instability and mass entrainment, 
the simulations we discuss here span a significant axial length and are considered completed only 
after much of the computed region reaches a quasi-steady state. This approach contrasts with 
that of numerical simulations with periodic boundary conditions along the jet axis, e.g., Bassett 
& Woodward (1995), Bodo et al. (1995), and Bodo et al. (1998), which focus on temporal growth. 
This paper is organized as follows: in §2, we explain the initialization of the 3D MHD simulations, 
and in §3, we discuss, in order, the growth of mass entrainment, the observational characteristics 
associated with mass entrainment in simulated (radio wavelength) intensity and polarization 
images, and the testing of the presence and measurable consequences of the KH instability. 

2. INITIALIZATION OF NUMERICAL SIMULATIONS 

These simulations were performed with the three-dimensional magnetohydrodynamics code 
ZEUS-3D, an Eulerian finite-difference code with the Consistent Method of Characteristics 
(CMoC) that solves the transverse momentum transport and magnetic induction equations 
simultaneously and in a planar split fashion (Clarke 1996). Interpolations are carried out by 
a second-order accurate monotonic upwinded time-centered scheme (van Leer 1977) and a von 
Neumann- Richtmyer artificial viscosity is used to stabilize shocks. The code has been thoroughly 
tested via MHD test suites as described by Stone et al. (1992) and Clarke (1996) to establish the 
reliability of the techniques. 

The useful computational domain for the three expanding jet simulations discussed here is 
a 3D Cartesian grid of 325 x 130 x 130 zones, which required roughly 105 Mwords of memory 
on the Cray C90 at the Pittsburgh Supercomputer Center (or PSC). The central thirty zones 
span the initial jet diameter, 2Rq, along the transverse Cartesian axes (the x-axis and y-axis). 
Outside the uniform grid zones, the grid zones are ratioed where each subsequent zone increases 
in size by a factor of 1.050. The 130 transverse zones span a total distance of 30i?o- Along the 
jet or z-axis, 225 uniform zones span the interval < z/R® < 30 and an additional 100 ratioed 
zones span an additional 30i?o with each subsequent zone increasing in size by a factor of 1.015. 
Outflow boundary conditions are used except where the jet enters the grid, where inflow boundary 
conditions are used. 

The three simulations are initialized at the inlet with parameter values similar to simulations 
in HCR and RHCJ. Specifically, simulation A here is equivalent to run C in RHCJ, simulation B is 
equivalent to run A, and simulation C is equivalent to run E. Simulations A and B contain jets that 
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are denser than the external medium and simulation C is of a light jet. Details of jet properties 
arc listed in Table [l|. While some simulations in RHCJ had a primarily toroidal magnetic field, 
all of the simulations here have a primarily axial magnetic field. Limiting this study to axial 
field jets was done to maximize the instability of the jet, since we anticipate expanding jets to be 
more stable and we know that a toroidal magnetic field can have an additional stabilizing effect 
(compared with a primarily axial magnetic field of the same strength, e.g., Appl & Camenzind 
1992; RHCJ). A smaller jet velocity, u, is used in order to reduce all of the Mach numbers (i.e., 
sonic, Alfvenic, and magnetosonic) because we wish the instability to grow significantly while 
the jet is within the computational domain and the growth length of the instability is roughly 
proportional to the product of jet radius and magnetosonic Mach number, with J\l ms = u/d ms . 
The magnetosonic speed is defined as o? ms = aj t + Vj, where the internal jet sound speed and 
the Alfven speed are defined by aj t = Tp/p and V\ = B 2 /Airp, respectively, and we assume the 
adiabatic gas law, P oc p r with T = 5/3. We have analyzed the simulations at a somewhat later 
time, at r = 45 dynamical times for simulations A and B compared with r = 36 in RHCJ and r 
= 16 for simulation C compared with r = 14 in RHCJ, where r = ta ex /RQ and a ex is the external 
sound speed. The jets in all of the expanding jet simulations pass through at least 3 flow-through 
times at z = 50i?o- Each simulation required about 40 CPU-hours on the Cray C90 at PSC. 

The jet is initialized across the computational domain with a uniform density, pjt, and initial 
radius, Rq. The magnetic field within the jet is initialized with a uniform axial component, B z , for 
r < Rq, where r is the radial position, and a toroidal magnetic component, -B^, with a radial profile 
that has the form B^r) = B^ k sin 2 [7r/(r)], where B vk is the maximum value for B^ and for r < r pk , 
f(r) = 0.5(r/r pk ) a , and for r pk < r < r max , f(r) = 1.0 - 0.5 [(1 - r/r max ) /(l - r pfc /r max )] 6 . In 
these simulations a = 1.106 and b = 0.885, which leads to a FWHM of the B^ radial profile equal 
to 0.5Ro- This radial profile for B^, which sets B^ to zero at v/Rq = 0.0 and for t/Rq > 0.9, is 
not physically motivated, but is mathematically well-behaved. The magnetic field is set to zero 
in the external medium. The jet thermal pressure is initially modified to satisfy an equation of 
hydromagnetic equilibrium, 



where the term on the right hand side describes the effects of magnetic tension. Note that this 
equilibrium condition is one of static balance, but does not necessarily lead to hydromagnetic 
dynamic equilibrium, which is determined by the condition V x (u x B) = 0. Owing to the small 
toroidal field in each of these simulations, this dynamic equilibrium is nearly satisfied and we have 
verified that no significant azimuthal velocity occurs within the jet near the inlet. The inlet values 
for B z and B^ k are given in footnotes to Table p] in the units used by ZEUS-3D, which sets the 
permeability of free space to 1. While the thermal pressure has been modified to satisfy the radial 
static equilibrium, the departure from a "top hat" pressure distribution is small. The typical 
plasma (3q = 8-Kpjt/B z is 200 in simulation A and 1.1 in simulations B and C. 
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In the simulations the external medium is isothermal and the external density, p ex (z), 
declines to produce a pressure gradient that is designed to lead to a constant expansion, 
Rjt(z) = (1 + z/120Rq)Rq, of a constant velocity adiabatic jet containing a uniform axial magnetic 
field and an internal toroidal magnetic field that provides some confinement, i.e., 



where g = R(z)/Ro, R{z) is the jet radius at axial position z, the constants C\ and C2 are based 
on the ratios at the inlet of jet thermal pressure to magnetic pressure from the axial component 
(in C\) or from the toroidal component (in C2), and C3 is a normalization constant equal to 1 + 
C\ - C2. We have confirmed that the initial expansion occurs quickly, within 5-7 dynamical times 
for the dense jet simulations and 1-2.5 dynamical times for the light jet simulation. 

A difference between these simulations and their cylindrical equivalents in RHCJ is a slightly 
reduced angular precessional frequency, co (see Table ||). We decrease u to obtain the same 
normalized precession frequency, uiR(z)/a ex , at some small (but nontrivial) distance from the 
inlet in each expanding jet as in the appropriate cylindrical equivalent. The spatial KH stability 
analysis of expanding jets suggests that at peak growth the normalized precession frequency is 
constant (Hardee 1986); our choice for the precession frequency, u, is below that for peak growth, 
but ujR{z)/a ex will increase with jet expansion. This avoids the unresponsiveness of a jet to too 
high a precessional frequency. The precession is in the counterclockwise sense at the inlet and 
should yield a clockwise spatial helical twist when viewed downstream from the origin/inlet. The 
simulations are analyzed at 6.5, 7.9, and 3.6 precessional periods for simulations A, B, and C, 
respectively. Another difference between the simulation sets is that the radial variation of the 
toroidal magnetic field component, Ba,, peaks at 0.5i?o i n the expanding jet simulations, while in 
all three cylindrically equivalent simulations in RHCJ this maximum in Ba, occurred at r = 0.8i?o- 
However, since the toroidal field is not dynamically important in these simulations, this difference 
is inconsequential. 



Since the code does not include a tracer of jet material, we use magnetic field and axial 
velocity criteria to estimate the mass associated with the jet plus any external medium that 
has been accelerated or mixed with the jet sufficiently to be called "entrained" (RHCJ). Such 
a velocity criterion has been used previously as an estimate for mass entrainment (Bassett &; 
Woodward 1995, Loken 1996). Since RHCJ found that the choice of the magnetic field or axial 
velocity criteria result in similar values for the estimated jet plus entrained mass, we have chosen 
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RESULTS 



3.1. Linear Mass Density as a Mass Entrainment Estimator 
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to use the magnetic field criterion alone in this paper. As in RHCJ, we define the linear mass 
density, a, at any point along the jet as 



where A is the cross-sectional area of the computational domain at z, and / is a switch set to 1 if 
the local magnetic field is above a threshold value and / = otherwise. We set the threshold to 
0.04B max , where B^ax = B^^z) 2 + B z (z) 2 . In a steady expanding flow, varies as RizY 1 and 
B z varies as R(z)~ 2 . While the choice of the 4% level is somewhat arbitrary, we believe that a 
threshold at this value is a useful demarkation between "mixed" (i.e., jet plus entrained material) 
and "unmixed" media. The unmixed external medium may contain some small magnetic field 
from numerical diffusion and should not be considered as entrained. In our previous simulations, 
this 4% B max level corresponded to external material accelerated up to ~ 15% of the jet speed in 
the dense jets and ~ 4% of the jet speed in the light jet. In a smooth flow of constant jet speed 
without entrainment, the linear mass density remains constant for both cylindrical and expanding 
jets. In what follows, <j{z) is the linear mass density for the axial position z normalized to the 
initial jet linear mass density. 

We display the spatial development of cr(z) for the expanding jet simulations and for the 
equivalent cylindrical jet simulations in Figure |l[ In addition, the bottom panel of Figure [l] shows 
the average axial velocity of the magnetized (with B > 0.04B max ) material, weighted by the linear 
mass density (i.e. (v z ) = J fpv z dxdy/ J fp dxdy), along z. From a comparison of the linear mass 
density in equivalent cylindrical and conical jet simulations, the jet plus entrained mass at the 
z = 60i?o boundary is reduced by roughly one-half in the dense jets and by at least two-thirds 
in the light jet (see Table ||). Calculations of the mean density of entrained material (see Figure 
^), which is defined as the mixed linear mass density less the initial jet linear mass density and 
corrected for different cross-sectional areas, shows that this mean density is similar to the initial 
external medium density. Thus, the primary reason for the reduction in entrained mass is the 
reduced external density in the expanding jet simulations, which is required by total pressure 
balance between the expanding jet and the isothermal external medium. For the parameters that 
we have simulated, the reduction of a in the expanding jets appears similar to (for the dense jets) 
or greater than (for the light jet) the reduction caused by changing the magnetic field configuration 
from primarily axial to primarily toroidal in the cylindrical jet simulations of RHCJ. 

In RHCJ, the linear mass density plots showed regions with different rates of increase of cr(z) 
with z, which we interpreted as different mass entrainment rates. Specifically, successive regions 
of slow, fast, and no mass entrainment (i.e., mass mixed with the fast moving jet as measured by 
our magnetic criterion) were found, and there was evidence that the three stages of growth were 
associated with linear, nonlinear and saturated growth of the KH instability. In both dense jets 
(simulations A and B), the rate of increase of linear mass density with axial position remains 
approximately constant for the entire range plotted. However, for z/Rq ;> 20 the fluctuation of a 




(3) 



-7- 



about this constant rate also increases with z, which is consistent with the notion of increasingly 
large "gulping" vortices, such as those discussed by DeYoung (1996), at the jet surface. This in 
turn is a signature of a developing KH instability. In addition, there are noticeable changes in 
dcr/dz near this position {z/Rq ~ 20) when the linear mass density is computed with a smaller 
threshold value (i.e., smaller than 4%). If the onset of the fluctuations in a in the dense expanding 
jet simulations is associated with the first transition point, then the first transition points of the 
dense jet simulations occur at roughly the same axial positions as that of the first transition points 
in their cylindrical equivalents (see Table ||). 

Of the dense expanding and cylindrical jet simulations, only the linear mass density in the 
cylindrical equivalent to the weak field jet (simulation A) could be interpreted as leveling off before 
the outer z boundary. From our analysis of the growth of magnetized mass in the simulations 
here and in RHCJ, we suspect that the linear mass density in the dense expanding jets would 
eventually saturate beyond the second transition (saturation) position in the cylindrical jets. If 
this is true, then the distance between the first and second transition points would be larger in the 
dense expanding jets than in their cylindrical equivalents. 

In the expanding light jet (simulation C), there are three regions of differing rates of increase 
of cr(z) with z. Estimated rates of increase are shown for the second (fast growth) stage in the 
light jet simulations in Figure |l[ Transition points between each region are listed in Table || for 
the expanding and cylindrical jets simulations. The first transition point is significantly closer to 
the inlet in the expanding jet simulation than in the cylindrical jet simulation. Since the second 
transition point is at roughly the same position in both the expanding and cylindrical jets, we see 
that the distance between the first and second transition points is larger in the expanding jet than 
in the cylindrical jet. As with the dense jets, the light expanding jet has a smaller spatial rate of 
increase of cr(z) with z in the fast growth region than the light cylindrical jet. 

We postulate that the KH instability is associated with the mass entrainment in these 
expanding jet simulations and compare the transition positions in Table ^ with estimates of the 
growth lengths of the KH instability determined from a linear analysis. Estimated growth lengths 
of the helical surface wave (n = 1, m = 0) for both the expanding and cylindrical simulations 
are listed in Table ||. All of these estimates are from the simple approximation for the growth 
length at the maximum growth rate, I* = /c^ _1 = — {2MJ^ S R)/ ln(Auj nm R/a ex ) (equation 5 in 
HCR), where Lo nm R/a ex = (n + 2m + l/2)7r/2 (equation 4a in HCR). Note that this method 
consistently underestimated the values found from more sophisticated root-finding techniques for 
the cylindrical jet simulations in RHCJ. Since the growth length depends strongly on the jet radius, 
for the expanding jet numerical integration was used to compute the number of effective e-folding 
lengths, N e , such that N e = J^ =Q ki(z)dz (equation [17] in Hardee 1986). From the smaller inlet 
magnetosonic Mach number, the expanding jets have a slightly smaller initial e-folding length 
than their cylindrical equivalents. The similar e-folding length for the expanding and cylindrical 
jets is consistent with the similar (or smaller in the light expanding jet) axial position of the 
first transition point. The increasing radius of an expanding jet leads to a larger distance for 5 
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e- foldings, which suggests that the expanding jet simulations should remain more stable farther 
down the grid. This greater stability is consistent with the larger distance spanned between first 
and second transition points in the expanding jets than in their cylindrical equivalents. 

The average axial velocity of the magnetized material decreases at the same rate with z in 
the dense expanding jet simulations (simulations A and B). This is similar to the constant rate of 
increase in the linear mass density, although the significant fluctuations in a are absent in (v z ). In 
the light expanding jet (simulation C), there are at least three regions where the average axial 
velocity behaves differently, and these regions are roughly coincident with the regions of differing 
mass entrainment rates. Specifically, there is an initial region out to z/Rq <^ 14 where there are 
very small fluctuations in the average axial velocity, which changes its smooth rate of decline at 
z/Rq ~ 8. Beyond this initial region with small fluctuations, there are significant fluctuations 
superimposed on a roughly constant decrease in the average axial velocity for 14 ^ z/Rq <J 35, 
and there are small fluctuations about an approximately constant (v z ) thereafter. The greater 
reduction in average velocity in the light jet simulation compared with the dense jet simulations 
is a result of mixing of the jet with the denser external medium. The region of approximately 
constant (v z ) in the light expanding jet suggests that any mixing has ceased and is consistent with 
the saturation of the KH instability. 

We calculated (but do not show) the transverse area covered by magnetized material with 
B > 0.04B max . In general, the area occupied by magnetized material is approximately the same 
for the dense expanding jets and for their cylindrical equivalents at any computed z. The area 
occupied by magnetized material in the light expanding jet simulation is slightly smaller than 
the area of its cylindrical equivalent. Therefore, the relative increase in the area occupied by 
magnetized material, i.e. (magnetized area)/(initial jet area), is smaller in all expanding jet 
simulations than in their cylindrical equivalents. This is consistent with a smaller spatial growth 
rate of the KH instability. In addition, we computed (but do not show) the linear momentum flux 
density, <jv z , of the mixed and unmixed regions. As in the trans- Alfvenic jets studied in Hardee 
&: Rosen (1999), the mixed material carries the bulk of the momentum flux. The maximum 
momentum flux density carried by the unmixed material in any computed plane of transverse 
zones is about 10% of the initial jet momentum flux density for the dense jets (simulations A and 
B) and about 5% for the light jet (simulation C). 



3.2. Comparison of Linear Mass Density and Simulated Radio Images 

In Figure || we display maps of simulated total radio intensity, which are line of sight 
integrations of Pth{B sin#) 3 / 2 , where 6 is the angle between the line of sight and the magnetic 
field (Clarke 1989). Such an approximation is necessary when relativistic particles are not tracked 
explicitly. The total intensity is integrated along the y-direction and is plotted on the xz-plane in 
the figure. In order to compare the morphology of the images, we have used different scales in 
each panel of Figure ||. In each panel, the range covers 3 orders of magnitude and the grayscale 
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maximum is 20% above the actual maximum intensity. The total intensity is overlayed by B-field 
polarization vectors. We note that B-vectors are aligned with filaments in the intensity image, 
as was also the case in the simulations in RHCJ. Of the three jets, the dense equipartition jet 
(simulation B) has the most obvious spine-sheath morphology in the intensity images and also 
maintains the spine farthest across the grid. The light jet simulation (simulation C) has a more 
uniform intensity for z/Rq ;> 40 than the dense expanding jet simulations, which is evidence for 
more uniform mixing and is consistent with a saturation of the KH instability. 

In the cylindrical jets studied in RHCJ, the slow, fast, and saturation stages of mass 
entrainment and the KH instability are associated with slow, fast, and zero expansion rates of the 
combined spine and sheath in the simulated intensity image. The transition points between the 
different expansion rates, as determined from Figure ||, are listed for both the expanding jets and 
their cylindrical equivalents in Table |E|. In the dense expanding jet simulations presented here 
(simulations A and B), a low intensity quickly expanding sheath appears at z/Rq ~ 20. This rapid 
expansion ceases for z/Rq ;> 50 in the simulation with the stronger field (simulation B), another 
indication of somewhat greater stability. The intensity image in simulation B appears to enter the 
third (saturation) stage although no equivalent stage for the magnetized mass (a) appears on the 
grid (see Fig. |l|). This suggests that this simulation continues to entrain mass for some distance 
beyond where the apparent expansion stops. This characteristic was seen in some of the cylindrical 
simulations of RHCJ, although this usually occurred for light cylindrical jets. In simulations A 
and B, the transition points of the total intensity expansion listed in Table || are similar to those 
listed for the linear mass density in Table |3|. We note that the fluctuations in the linear mass 
density may be associated with growing vortices at the jet/external medium interface, and there is 
evidence for such vortices in the wispy edges of the total intensity images in simulations A and B. 

There are three stages of expansion in the total intensity image from the light jet simulation 
C; particularly notable is the nearly constant width for z/Rq ^ 45. Again we see a correspondence 
between the transition positions determined from the linear mass density and from the total 
intensity width. As with the transition points demarking different mass entrainment rates, the 
distance between the two transition points determined from the intensity images for the expanding 
jets is larger than the distance between the transition points for the cylindrical jets. This result 
is consistent with the initially shorter e-folding lengths but longer distances to 5 e-foldings on the 
expanding jet when compared to the cylindrical jet (see Table |j). Thus, jet expansion has served 
to partially stabilize the lower magnetosonic Mach number expanding jet relative to a cylindrical 
jet and the total intensity images show observational consequences of this stabilization. 

3.3. Jet Structure Related to Mass Entrainment and KH Instability 

In this subsection, we compare wavelengths and wave speeds measured in the simulations 
with those expected from a linear analysis of KH induced modes. Subsequently, we examine the 
effect of the KH instability on the fluting of the jet surface via axial velocity cross-sections and on 
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internal structure via axial cuts of each velocity component. 

3.3.1. Wavelengths and Wave Speed Estimates 

From the oscillations observed in the intensity images (Fig. |3|), we estimate the wavelength of 
the surface helical mode, A/j, to be: in simulation A, 15i?o; in B, 13-Roi an d in C, 18i?o- The error 
in this estimate is roughly 10%. We measured these wavelengths in a region centered on z/Rq = 
25 in simulations A and B, and z/Rq = 20 in simulation C; recall that the wavelength at maximum 
growth should be proportional to the jet radius. The expanding jets have a smaller precessional 
frequency, a lower jet velocity and a different wave speed than their cylindrical equivalents. These 
factors combine to give similar wavelengths in the expanding jets and their cylindrical equivalents 
(see Table ||) . 

If the observed wavelength of a mode is much larger than the wavelength associated with peak 
growth, A/i > 10A*, then the wave speed should be roughly rju/(l + rj). However, if A* < A^ < 3A*, 
the wave speed approximations for peak growth are more appropriate, and the wave speed should 
be ~ r/ 1 / 2 u/(l + 77 1 / 2 ) (e.g., Hardee 1987, and HCR). From the approximation (4b) in HCR for 
the wavelength at maximum growth, we estimate in simulations A and B at z/Rq = 30 that 
A* = 10i?o and in simulation C at the same position A* = 16i?o- Thus, the measured wavelengths 
indicate that the wave speed in these simulations should be ~ T] l l 2 u/(1 + in 1 / 2 ). Assuming that 
the initial values for the external density and the jet densities follow from the expected expansion, 
this approximation yields the wave speeds given in Table § These wave speeds are consistent with 
wave speeds estimated from the product of the observed helical mode wavelength and the driving 
frequency, v = uj/2-k (also in Table |6|) . 

In order to calculate the wave speed directly, we measure the movement over time of 
maxima in axial profiles of the simulated total intensity images (Figure ||). The linear analysis 
approximation for wave speed suggests that the wave speed should increase only moderately with z 
in the dense jet simulations. On the other hand, the wave speed should increase more dramatically 
in the light jet simulation. Specifically, in the weak magnetic field dense jet (simulation A), the 
wave speed increases from ~ 2.3a ex at z/Rq ~ 15 to ~ 2.7a ex in the 25 <^ z/Rq ^ 40 interval. In 
the more strongly magnetized dense jet (simulation B), the wave speed increases from ~ 2.7a ex at 
z/Rq ~ 10 to ~ 2.9a ex for 30 <; z/Rq <J 50. Thus, the wave speed in the dense jets increases by 
roughly 10-20% over this range of axial positions. In contrast, in the light jet simulation the wave 
speed increases from ~ 3.0a ex at z/Rq ~ 15 to ~ 7.6a ex at z/Rq ~ 30. The large wave speed, 
which occurs where the KH instability has reached nonlinear proportions, does not increase much 
beyond this axial position. It is likely that this constant wave speed is related to saturation of 
the KH instability (as inferred from the saturation in the linear mass density and the constant 
apparent width of the intensity image) for z/Rq ^ 40. 

The measured wavelengths and wave speeds, except where the amplitude of the KH instability 
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has grown beyond the linear approximation, are consistent with those estimated from the linear 
analysis of the KH instability. Thus, we conclude that the KH instability is responsible for the 
growth of mass entrainment in these simulations. 

3.3.2. Cross-Sections of Axial Velocity 

The growth of the KH instability and its effect on the jet is dramatically demonstrated in the 
grayscale cross-sections of axial velocity from simulations A, B, and C in Figure |5|. As observed 
in magnetic pressure cross-sections in previous simulations (RHCJ), there are many corrugations 
in the fast moving surface at z/Rq ~ 15 in all three simulations. That the many corrugations are 
roughly evenly spaced circumferentially on the jet is evidence for the surface waves of high order 
fluting modes. In all the simulations here and in simulations in RHCJ, the higher order modes that 
dominate close to the inlet are overwhelmed by the larger distortion amplitudes accompanying the 
slower growing, lower order modes farther down the jet. In Figure [| the helical mode appears as 
the clockwise motion of the jet center about the initial jet axis as one moves down the jet. This 
clockwise rotation about the jet axis provides an additional estimate for the wavelength of the 
helical mode. In the light jet simulation, Xh/Ro ~ 35 between z/Rq = 25, where much of the jet is 
to the left of the original jet axis, and z/Rq = 40, where the jet has been displaced to the upper 
right. This wavelength is consistent with the very large wave speed (v w /a ex = 7.6) measured at 
large z/Rq. 

There is also a correspondence between features in the intensity images and the axial velocity 
cross-sections. Recall that the simulated intensity images are integrated along the y axis and 
the observer in Figure |3] is to the right in Figure ^. For the light expanding jet simulation, 
a spine-sheath morphology is more evident in the cross-sections than in the intensity images. 
However, far from the inlet the area of the high v z material is smaller than for the dense expanding 
jets. This is another indication of greater mixing occurring in the light jet simulation. 

3.3.3. ID Axial Cuts of Velocity 

The effect of the KH instability on the internal structure of the jets is demonstrated by the 
one-dimensional cuts of velocity components along the jet axis (x,y = 0) shown in Figure ^. In all 
three simulations, the variation in the axial velocity component is initially small, <^ 2%. At some 
point the variation grows to significant amplitude relatively abruptly. For the dense jet simulations 
this occurs at z/Rq ~ 22, and for the light jet simulation this occurs at z/Rq ~ 15. These positions 
are close to the first transition point (in either magnetized mass or width of the simulated intensity 
image) in all three jets. Thus, changes in the nature of the mass entrainment rate associated with 
the growth of the KH instability are able to affect the jet dynamics significantly. 

Continued growth of the KH instability as illustrated by axial profiles of v z differs in the dense 
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jet simulations beyond the first transition point. Specifically, there are significant oscillations of v z 
in simulation A for 20 <J z/Rq <J 42, but oscillations in simulation B remain relatively small out 
to z/Rq ~ 47. The smaller variations in the more strongly magnetized jet imply greater stability. 
In the light jet simulation, there is an additional, extremely abrupt (0.75u) decrease in v z at 
z/Rq ~ 31. Examination of the axial velocity cross-sections in Figure || reveals that this extremely 
abrupt drop in axial velocity in the light jet simulation is created as the small high velocity spine 
is displaced off the axis by the helical mode of the KH instability, and is not associated with the 
onset of the plateau in the jet plus entrained mass. The second transition position in the light jet 
does coincide with the reappearance of the v z maxima in Figure ||. 

In the dense jet simulations the transverse velocity components show a noticeable relative 
phase shift in the region 10 <^ z/Rq <^ 30. We consider the maximum in v y at z/Rq = 10 in 
simulation A to be a precursor of, and therefore connected to, the maxima in v x at z/Rq = 21 
and in v y at z/Rq ~ 23. The phase shift is more easily seen in the more strongly magnetized jet 
(simulation B), with extrema in v x at z/Rq = 12, 18, and 23 repeated in v y at z/Rq = 13, 21, 
and 25. In the light jet simulation, the evidence for a relative phase shift between the transverse 
velocity components is less obvious, but perhaps the maximum in v x at z/Rq = 7 contributes 
to the maximum in v y at z/Rq = 10. In all of the simulations, the extent of the region where 
the phase shift is noticeable is related to the growth of the helical mode of the KH instability: 
the region begins where the fluctuations in transverse velocity reach a significant amplitude and 
ends where the displacement from the KH instability has become significant. Where noticeable, 
v y is shifted farther down the jet than v x . Such a sequence is in the sense of a clockwise spatial 
variation when viewed toward +z from the inlet and consistent with the counterclockwise temporal 
precession at the inlet described in §2. This phase shift is additional evidence for a significant 
presence of the helical KH mode. 

4. CONCLUSIONS 

The magnetized, or jet plus entrained, mass in these expanding jet simulations increases 
more slowly than in their cylindrical equivalents. This reduction in the mass entrainment is 
a direct consequence of the reduced external mass density, which is required by an expanding 
pressure-matched jet. For the parameters that we have simulated, the reduction of magnetized 
mass in the expanding jet simulations appears similar to (for the dense jets) or greater than (for 
the light jet) the reduction caused by changing the magnetic field configuration from primarily 
axial to primarily toroidal in the cylindrical jet simulations of RHCJ. Thus, a conical jet appears 
more stable than its cylindrical equivalent, and should be considered, along with relativistic jets 
(Hardee et al. 1998; and Hughes, Miller, & Duncan 1999), as another stabilizing influence in the 
propagation of astrophysical jets. 

As in the cylindrical jet simulations, the spatial development of the magnetized mass passes 
through different stages, in which the spatial rate of increase of magnetized mass or the smoothness 
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of this increase is different. We see evidence for a first transition, between a region of a nearly 
constant increase in magnetized mass and a more varied increase in magnetized mass in the 
dense jets. In the light jet simulation this first transition occurs between a region with a slow 
increase in magnetized mass to one with a faster increase. Also, in the light jet simulation we 
see a second transition, between the region of fast increase and a region where the magnetized 
mass remains relatively constant. In the expanding jet simulations the distance between first and 
second transition points is similar to or larger than this distance in the cylindrical jet equivalent 
simulations. The observed behavior is consistent with slower growth of the KH instability on 
expanding jets relative to cylindrical jets. Thus, we confirm the theoretical prediction that the 
KH instability grows on these expanding jets at a reduced rate relative to an equivalent cylindrical 
jet. Previous work (RHCJ) associated the different stages of spatial development with the linear, 
nonlinear, and saturated stages of the KH instability, and this appears to be the case in these 
expanding jet simulations. The simulations also contain wavelengths and wave speeds that are 
consistent with theoretical estimates of the wavelengths and wave speeds appropriate to the 
surface helical mode wave triggered at the precession frequency. 

Simulated intensity images also reveal three stages in the expansion of the mixed jet and 
entrained material. These stages are similar (slow, fast, and zero expansion) and roughly coincident 
with the three stages of mass entrainment. However, in one case the jet continues to entrain mass 
while maintaining a constant jet width in the intensity image. The reduced KH instability growth 
rate in expanding jets does result in a lengthening span between transition points as determined 
from the intensity image relative to transition points for cylindrical equivalent jets. 

Our present study confirms results from RHCJ. The simulated total radio intensity images 
and the axial velocity transverse cross-sections show that dense jets are able to maintain a high 
velocity spine as part of a spine/sheath structure. We also see an alignment of polarization vectors 
with filaments in the simulated intensity images. In addition, the transverse cross-sections show 
a progression from high order, fast growing, small amplitude to low order, slower growing, large 
amplitude KH modes. Axial cuts of each velocity component show that the internal structure 
exhibits a transition in the growth of the instability in roughly the same position as the first 
transition point in the mass entrainment rate or in the width of the intensity images. Beyond this 
transition point, especially in the light jet simulation, the displacement of the jet from the initial 
axis complicates the profile. 

A. Rosen and P. Hardee acknowledge support from the National Science Foundation through 
grant AST-9802955 to the University of Alabama. The authors would also like to acknowledge 
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and the storage facilities associated with the Cray T90 at the San Diego Computing Center. 
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van Leer, B. 1977, J. Comput. Phys., 23, 276 

Fig. 1. — The normalized linear mass density a (top and middle panels) and average axial velocity 
(v z ) (bottom panel). Also shown is a from the cylindrical equivalent simulations in RHCJ. In order 
to show the extent of the second (fast growing) stage in the light jet simulations, estimates of the 
mass entrainment rates have been added in this stage to the middle panel. 

Fig. 2. — The average density of the entrained mass, p en , in the expanding jet simulations 
normalized by the entrained mass density in each equivalent cylindrical jet simulation. The solid 
line is the initial profile of the external density, the dotted line is the calculated mean density of 
magnetized mass less the jet mass. 

Fig. 3. — Integrated simulated intensity images overlayed with B-field polarization vectors. The 
total intensity is shown in grayscale, and covers three orders of magnitude with a maximum set 
to roughly 1.20 times the actual maximum intensity in the image. In addition, the polarization 
vectors have a length proportional to the fractional polarization and are shown only where the total 
intensity is above 0.001 of the scaled maximum intensity in the image. 

Fig. 4. — Evolution of total intensity profiles on the z-axis. The simulated total intensity, which is 
given in a quasi-logarithmic scale, is displayed for three different times for each simulation. 

Fig. 5. — Grayscale cross-sections of axial velocity in the xy-plane for simulation A (top), B 
(middle), and C (bottom). Here, the +x axis is towards the top of each panel and +y to the right 
(+z is into the page), and each axis spans the region -4 < x/Rq, y/Ro < 4. Darker shading indicates 
larger v z , white indicates v z < 0. The number in the upper left of each panel indicates the axial 
position, z/Rq. 

Fig. 6. — One dimensional cut of each velocity component along the z-axis. The panels on the left 
show the axial profile of v z and the panels on the right show v x and v y , all in units of the external 
medium sound speed, a ex . Note the different ranges of velocity on each panel, although this range 
is approximately 0.25w in all of the panels. 
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Table 1. Inlet Jet Parameters 



Sim. 


(vT 


(B) b 


(Pjt/PexY 


(a jt /a ex ) c 


(V A /a ex ) d 




) c ' d M ex 


(M jt f 


(M A ) d 




A 


4.00 


W 


0.995 


0.50 


0.04 


0.50 


3.49 


7.00 


93.1 


6.98 


B 


4.00 


E 


0.54 


0.37 


0.38 


0.52 


3.67 


10.00 


9.80 


7.00 


C 


0.25 


E 


0.54 


1.47 


1.50 


2.10 


11.17 


7.60 


7.45 


5.32 



& r\ = pjt/pex, uniform jet density across the jet. 

b Descriptions of the magnetic field: E=Equipartition (B z = 0.75, B^ k = 0.26, in ZEUS-3D units), 
W=Weak (B z = 0.075, Bj* = 0.035, in ZEUS-3D units). 

c Jet thermal pressure used is before slight modification (based on Equation [1]) to a uniform profile 
across jet. 

d Based on axial magnetic field only, valid for r = and 0.9 < t/Rq < 1.0. 



Table 2. Precessional Data 



Simulation u;R/a ex (ojR/a ex ) cy i a t p (v t /a ex ) b {v t /a™ Q ) h 



A 


0.9 


1.0 


6.98 


0.019 


0.039 


B 


1.1 


1.2 


5.71 


0.020 


0.039 


C 


1.4 


1.6 


4.49 


0.062 


0.030 



frequencies used in previous cylindrical jet simulations (RHCJ). 

b Only the transverse velocity related to the precession, does not 
include velocity associated with the jet expansion. 
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Table 3. Magnetized Mass Related Quantities 





a(z/R = 60) 


Transition Points 


Simulation 


Conical Cylindrical 3, 


Conical Cyl. a 


A 


1.75-2.0 2.8/3.6 b 


20-25 c />60 25-27/48 


B 


1.5 2.8 


20-25 c />60 25/>60 


C 


7-8 25 


10/42 18/40 



Cylindrical jet data from RHCJ. 

b The two values listed here depend on whether the jet reached 
saturation at z/R = 60 (2.8) or not (3.6). Both interpretations 
are possible. 

c Based on position where ever increasing fluctuations in a begin. 
Note that similar fluctuations occur at z/Rq ~ 20 as well in 
cylindrical jet simulations. 



Table 4. Effective Growth Lengths 



Expanding Jet Cylindrical Jet 
Simulation N e = 1 N e = 5 N e = 1 N e = 5 a 



A 


6.5 


39.5 


6.9 


34.5 


B 


6.5 


40.8 


6.9 


34.5 


C 


4.9 


28.9 


5.3 


26.5 



a Values here are 5 times those in previous column. 
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Table 5. Data from Simulated Intensity Maps 





Transition Points 






Simulation 


Conical 


Cyl. a 


Conical 


Cyl. a 


A 


20/>60 


25/45 


15 


13 


B 


20/50 


30/50 


13 


13 


C 


15/45 


-20/40 


18 


20 



Cylindrical data taken from equivalent simulations 
in RHCJ. 



Table 6. Estimated Wave Speeds 



Simulation 


v w {z/R = 0)/a ex 


v w (z/R = 20)/a ex 


\ h uj/(2ira ex ) 


A 


2.3 


2.4 


2.2 


B 


2.5 


2.6 


2.3 


C 


3.7 


4.0 


4.0 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 



Simulation A (WD) 
Simulation B (ED) 
Cylindrical equivalent to A 
Cylindrical equivalent to B 





30 

20 

10 




10 
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Simulation C (EL) 
Cylindrical equivalent to C 



First transistion for Sim C 
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